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Abstract 

Recent work on the theme of power corrections in perturbative QCD 
is briefly reviewed, with an emphasis on event shapes in e'’'e“ annihi¬ 
lation. The factorization of soft gluon effects is the main tool: it leads 
to resummation, and thus highlights the limitations of perturbation 
theory, pointing to nonperturbative corrections whose size can be es¬ 
timated. Power corrections can be resummed into shape functions, 
for which QCD-based models are available. Theoretical progress is 
closing in on the nonperturbative frontier. 
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1 Introduction 


It has been known since the early days of QCD [Q] that the perturbative 
expansion for IR safe observables is at best asymptotic, and in fact not even 
Borel summable. This is good news: it would be very surprising if perturba¬ 
tion theory alone could give a well-dehned answer in a theory such as QCD, 
where nonperturbative phenomena govern the spectrum of physical states; 
an ambiguity in the perturbative answer tells us that conhnement physics 
must, at some level, be present and relevant. Under the mild assumption 
that QCD is a consistent theory, we can do more: gauging the size of the 
uncertainty in the perturbative answer, we can estimate the impact of the 
dominant nonperturbative effects. 

The basic tools for this analysis are the factorization and the resummation 
of soft gluon effects. Soft gluon emission has a universal character, and 
factorizes from the hard part of scattering amplitudes and cross sections; 
factorization in turn implies that the dominant (logarithmic) contribution 
of soft gluons can be computed to all orders. These computations display 
explicitly the asymptotic nature of the perturbative series and can be used 
to study power-suppressed corrections. 

Studies of power corrections with perturbative methods have generated a 
vast literature in the past several years p[, and have reached an impressive 
degree of phenomenological success. In this short review I will focus mostly 
on theoretical results obtained for event shapes in e’''e“ annihilation, trying 
to put them in the perspective of the most recent developments. 

2 Perturbative windows on power corrections 

QCD resummations typically yield expressions of the general form 

Uq^) = ^{kra.(e) ■ ( 1 ) 

Such expressions are ill-dehned in perturbative QCD, because of the Landau 
pole in the running coupling at = A^. Expanding in powers of ctsiQ"^), one 
hnds that the coefficients of the expansion grow factorially. The size of the 
ambiguity in the result can be gauged by taking the residue of the pole: one 
hnds that it is suppressed by a power of the hard scale, Sfa/fa oc (A^/g^)“. 
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A classic example [Q is the resummation of fermion bubble insertions into 
a gluon line. Let n(fc^) denote the fermionic part of the gluon vacuum polar¬ 
ization, and let a(x) be a generic (possibly weighted) cross section depending 
on some kinematic variable x. Summing over insertions of If in the single 
gluon contribution to a gives an expression of the form 

^ J 'W |i + n(p)|2'^ \ V j ’ ^ 

where a is the virtuality distribution of the emitted gluon. If the rif de¬ 
pendence of the answer is reinterpreted as the abelian contribution to the 
running of the coupling, Eq. displays the Landau pole in the integration 
over gluon virtuality. 

The resummation of multiple soft gluon emission near the boundary of 
phase space leads to similar conclusions. Consider, for example, the Laplace 
transform of the thrust distribution |Q]. It can be shown to exponentiate in 
the form 

= (3) 

Jo dt ’ ^ ^ 

where the function S contains all the singular logarithmic dependence on 

t = 1 — T, and can be written (to NLL accuracy) as 




d/c I 


( 4 ) 


In this case the Landau pole obstructs the integration over transverse mo¬ 
mentum. 

A variety of regularizations for the Landau singularity have been pro¬ 
posed, (IR cutoff IQ , principal value prescription Q , regular IR continuation 
of the running coupling Q, dimensional regularization Q, choice of contour 
in the inverse Laplace transform Q). As far as the parametric size of the 
power corrections (he. the power of q with which they scale) is concerned, 
all these prescriptions (sometimes after some debate) lead to the same con¬ 
clusion. They may differ when more detailed prediction are attempted. 

A benchmark to verify the reliability of the methods based on different 
resummations is the application to observables for which nonperturbative 
information is available through other means, typically because of the appli¬ 
cability of the OPE. Along these lines resummation and OPE methods have 
been shown to give compatible results for the size of the correction in all 
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tested cases (ranging from the beautiful results of David on the nonlinear 
a model in d = 2, to the classic analysis of the total annihilation cross sec¬ 
tion by Mueller |^, to the more recent results on deep inelastic structure 
functions H). Making more detailed predictions, going beyond the iden- 
tihcation of the leading power correction, it is necessary to rely upon the 
assumption of ultraviolet dominance To introduce the idea, consider the 


OPE for an inclusive cross section, say a DIS structure function in Mellin 
space 




VP(iV,F/,9"){Or>(iV,F/)) + . 


q- 


(5) 


Here the factorization scale /i/ is an IR cutoff for the coefficient functions 
Oj, but it is an ultraviolet cutoff for the operator matrix elements. Physical 
quantities must not depend on /xj, and in fact the cancellation of the logarith¬ 
mic dependence on /ij, within a given twist, is enforced by Altarelli-Parisi 
equations. One observes, however, that the coefficient functions Ci have an 
ambiguous power dependence on /x/, of IR origin, due to the divergence of 
their perturbative expansion. At the same time, higher twist operator matrix 
elements have power-like UV divergences, so that they mix under renormal¬ 
ization with lower twist operators. These two facts make the twist separation 
ambiguous, and the ambiguity can only be disentangled once the same reg¬ 
ularization prescription is chosen for both Q and (Oj). What renormalon 
models and resummations are trying to do is to predict the exact dependence 
of higher twist operator matrix elements on kinematical variables (such as N 
above), based on their renormalization properties, i.e. on their UV behav¬ 
ior. These prediction will be correct to the extent that the matrix elements 
are “UV dominated”. Evidence supporting the idea of UV dominance has 
recently been presented in the case of DIS structure functions, in the 
elastic limit x ^ 1. In that limit matrix elements of operators of twist 4 and 
6 are dominated by parton conhgurations mimicking twist 2, so one can ex¬ 
plicitly verify that the N dependence of the perturbative ambiguity at twist 
2 (and large N) is in fact cancelled by the corresponding ambiguity of the 
higher-twist operator matrix elements. 

To go beyond the OPE, it is necessary to either rely upon renormalon 
models, or to take the viewpoint suggested by the factorization properties 
underlying soft gluon resummations. In the following, I will use factoriza¬ 
tion to dehne a way to parametrize power corrections to event shapes (the 
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shape function), following [rB|; starting from this general parametrization it 
is possible to recover the results of renormalon calculus, which in turn can 
be viewed as a QCD-motivated model for the shape function. 


3 Two—jet limit and shape functions 

Generic IR safe observables in production processes cannot be described with 
the OPE: they are weighted cross sections. Specihcally, consider the distri¬ 
bution of an event shape e, chosen so that the two-jet limit correspond to 
e —*• 0 (for example, e = 1 — T, G, pj,...). If the explicit expression for e 
in terms of the momenta of hnal state particles is e = E^ipi, ■ ■ ■ ,Pm), one 
dehnes 

? = / dLlPSn, \Mmf S{e- Em{Pi, ■ ■ ■,Pm)) , (6) 

ae Zq ^ J 

as well as the ‘radiation function’, R{e) = Jq de'da/de'. 

Soft gluon effects dominate the cross section in the two-jet region. Final 
state gluons are forced to be either soft, or collinear to the two back-to-back 
jets; as a consequence, the distribution factorizes under a Laplace (or Mellin) 
transform, schematically as (T(e, = Ji{q^e)*J 2 {q^e)*S{q^e^)*H{q^). There 

are, in fact, two relevant mass scales in the e —> 0 limit: the squared invariant 
mass of the jets, which vanishes as g^e, and the (squared) total energy carried 
by soft gluons, which vanishes as g^e^. The factorization is valid, within 
perturbation theory and for e —> 0, in the range q^ » q^e » q^e^ >> A^. 
If we assume that we can rely upon the same factorization in the extended 
range g^ >> g^e >> g^e^ ~ A^, we can use the perturbative results to 
parametrize power corrections of the type (A/(ge))^. 

Physically, the picture emerging from factorization is equivalent to the 
phenomenological ‘tube model’ |T^: neglecting inverse powers of g^e, one 
observes that jet masses are insensitive to the transverse components of soft 
gluon momenta, and are linearly shifted by an amount proportional to the 
total light-cone momentum carried by soft gluons. 


Mt 


NP 


= M? 


PT 


+ e^g , 


(7) 


where i labels the two jets and Cj is the sum of soft gluon momentum com¬ 
ponents along the i-th jet direction. At least for event shapes vanishing with 
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Mf, the effect of soft gluons on the perturbative distribution will take the 
form of a convolution in the light-cone components e^. Introducing an IR 
cutoff fi, one writes 

R{e) = j deide2f{ei,e2)RpTie) , ( 8 ) 

-'O Mf^Mf+etq 

where /(ei,e 2 ) is the most general form of the announced shape function. 
Note that it does not depend on the hard scale q. For shape variables such 
as C or t, which depend only on the sum of two jet masses in the two-jet 
limit, one can simplify the convolution by introducing 

/(e) = / (ieide 2 /(ei, 62 ) 5(6 - ei - 62 ) . (9) 

JO 

Factorization further provides an explicit, if formal, operator expression for 
the shape function /. In the limit we are considering, in fact, soft gluon 
emission is well approximated by replacing the two back-to-back jets by 
eikonal lines. One can then treat the joint distribution of 61,2 as a generic 
event shape, according to Eq. (H), but replacing the matrix element Adm 
with its eikonal counterpart. Essentially 

/('■.«=) = 4 E / ‘'LIPS" \MtfS («, - kr“] S (e, - kf') , (10) 

m '' 


where dLIPS^'^ is the m-particle phase space with hxed eikonal lines, and 
is the total soft gluon light-cone momentum flowing into hemisphere i. 
Such eikonal cross sections are matrix elements of Wilson lines, so the shape 
function can ultimately be expressed in terms of correlators of the energy- 
momentum tensor at large distances, in the presence of eikonal sources. 

To compare the results of this approach with renormalon calculations, it is 
useful to see how the shape function emerges from a resummed perturbative 
calculation. Consider the expression for thrust, Eq. (^. There, small t is 
associated with large u, and the dominant contributions to the integral arise 
from the region a \/v t. The leading power corrections that we are 
trying to resum are then related to the IR behavior of the coupling through 
the cusp anomalous dimension F, and one can isolate them by interchanging 
the order of integration, and introducing a transverse momentum cutoff. One 
dehnes 


SNp{u/q,fj) 


/■m" dkl 

Jo 


f /; 2 n \ da 
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= W). 


n=l 


n\ 


( 11 ) 


where terms suppressed by powers of v jq^ have been neglected and 

dklk’l^-T (a.(kl)) ( 12 ) 

are an inhnite set of dimensionful nonperturbative parameters, controlling 
power corrections of decreasing size, but always of the desired form 
They can be organized into a shape function via a Laplace transform 

exp (-S'Arp(z//g,/r)) = / ft{e, fi) . (13) 

Jo 

Within this perturbative framework, it is easy to make contact with the 
dispersive approach P] , whose results are recovered at leading power. Specif¬ 
ically, 


the dispersive approach predicts moments of event shapes in terms of a 
single nonperturbative parameter aQ. This parameter is equivalent (in 
fact, in a suitable factorization scheme, proportional) to Ai. From the 
point of view of the shape function, one can show that the parameters 
A„ with n > 1 contribute to power corrections to moments only at the 
level of subleading corrections, 1/q^ with p > 1. Power corrections for 
the moments of the distribution are then, in fact, well approximated 
by retaining only Ai: it is only in the deep nonperturbative regime 
e ~ A/g that subleading power corrections play a signihcant role. It is 
interesting to note [|I^ that within the dispersive approach central mo¬ 
ments (such as — (t)^) receive no leading power corrections, 

whereas these are nonvanishing in shape function hts. 


In the dispersive approach one needs to include a correction accounting 
for the correlations between emissions in different hemispheres O, 12 , 


h7|, ITSl . The shape function accounts for this effect by a lack of factor- 
izability: if one writes /(ei, 62 ) = g{^i)g{^ 2 ) + ^/(ci, ^ 2 ), the 6f term is 
responsible for correlations. Clearly, renormalon calculations provide 
a model estimate of this effect, whereas from the point of view of the 
shape function this will emerge as a result of hts to data. 
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Tube and renormalon models predict that the leading power correction 
shifts event shape distributions by an amount proportional to the en¬ 
ergy carried by soft gluons. This prediction is recovered by the shape 
function, which introduces additional smearing due to subleading power 
corrections, as easily seen by using eq. |T^, and as described in detail in 
Ref. ^ 
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4 Dressed gluon exponentiation 

It is clearly desirable to merge the information that can be extracted from 
soft gluon resummation (which expresses the dynamical effects of pertur¬ 
bative multiple soft gluon emission) with the results of the dispersive ap¬ 
proach (which models the effects of nonperturbative single gluon emission). 
This merging is achieved by the method of dressed gluon exponentiation 
(DGE). 

To summarize the main features of the method, pick as a shape variable 
a jet mass p. Schematically, one proceeds as follows . 


The hrst step is the computation of the characteristic function of the 
dispersive method (that is, the differential cross section for the emission 
of a gluon with virtuality fc^) for the observable at hand, in the Sudakov 
limit. In this limit, one retains only terms that are singular as p —> 0 
when k’^/{pq^) is kept hnite, as these are the only terms contributing 
to the Sudakov logarithms in the two-jet limit. For the jet mass one 
hnds, setting ^ 


r{p,i) 


log 


2 

P 


l_e 

p^ 


(14) 


Within the dispersive approach, the characteristic function must be 
integrated over virtuality with a weight given by the running coupling. 
Here one encounters the Landau pole. One possibility to dehne the 
integration is to turn to a Borel representation of the coupling |]2T |. 
In the present case, to achieve NLL accuracy it is crucial to dehne 
the coupling so that it satishes Altarelli- Parisi evolution at NLO, at 
least for the singular contributions near threshold. This is the “gluon 
bremsstrahlung” coupling 0- Performing the integration over gluon 
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virtuality, the result for the distribution is of the form 


^ = duB{u, p) exp ^^^^Ab{u) . (15) 

dp /po JO ^ J TTU 

Here Ab{u) is the Borel transform of the chosen coupling, the factor 
sin7rM/(7rM) arises from having taken the timelike discontinuity of the 
coupling, as prescribed by the dispersive method, and B (u, p) is the 
Borel function obtained by integrating Eq. (|I^ ) over virtuality with 
the appropriate weight. 


The key step to tie this single gluon result to resummation is to use 
the dressed gluon distribution as kernel of exponentiation. One writes 
the Laplace transform of the jet mass distribution as 

In ^ ^ - l) dp , (16) 

and substitutes in this expression the single dressed gluon result, given 
by Eq. (0). 

At this point, the Borel representation of the exponent suggests a pat¬ 
tern of exponentiated power corrections. The integral transform asso¬ 
ciated with the exponentiation turns the simple pole structure of the 
one-loop Borel function B (u, p) into a much richer one for the trans¬ 
formed function B{u, v). One can use this structure to model the pat¬ 
tern of power corrections. The radiation function can be reconstructed 
by inverting the Laplace transform, using 




dv 


■ exp 


/c 2nii> 

with a nonperturbative contribution 


z/p + lnJ^^(z/, g) + lnJ^^(i/, q) 


(17) 


1,19) 


C50 -1 

EAn- 

n=l 


n\ 



(18) 


In this approach, A„ are still free parameters, however the structure of 
the Borel transform may £x relations between them, or other general 
properties. For example, in the case of the jet mass p one is led to 
conclude that \ 2 k = 0. 






• These results can again be interpreted in terms of a shape function, 
writing g) as in Eq. (|l^). Now however the structure of the 

moments is derived from the Borel representation. 


Dressed gluon exponentiation provides a useful tool to merge the predictions 
of soft gluon resummation with a parametrization of power corrections dic¬ 
tated by renormalons. The results of resummation are reproduced to NLL 
accuracy, but further predictions are made: in fact, the size of all subleading 
perturbative logarithms is predicted in the large Uf limit, and found to grow 
factorially. This fact, perhaps not unexpected, in view of the renormalon 
singularities of the original resummed expression, sets limits on the validity 
of resummations performed in terms of successive towers of logarithms, and 
provides an independent tool to estimate the range in which they are appli¬ 
cable. As far as power corrections are concerned, DGE shows that in the 
two-jet limit at least a class of them exponentiates together with Sudakov 
logarithms. Further, DGE provides a dehnite regularization prescription to 
handle resummed perturbative expressions at power accuracy, which may be 
useful to check the assumption of ultraviolet dominance whenever the explicit 
operator form of power suppressed contributions is known, as was done for 
DIS in Ref. [M. Finally, as announced, renormalons make specihc predic¬ 


tions on the structure of the shape function for a given observable, which in 
principle can be tested against data. 


5 Phenomenology 


The shape function approach and DGE were applied to selected event shapes 
0 d [^]. For lack of space, I must refer the reader to the original 
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papers for details, however some observations are in order. 

First of all, as pointed out above, DGE puts constraints on the form 
of the shape function, which may depend on how renormalon calculus is 
implemented. For example, Korchemsky and Tafat, in Ref. |T^, include in 
their £t the effect of correlations between hemispheres by picking a gaussian 
ansatz for the shape function with a non-factorizable term; this term is 
found to be important for the fit. Gardi and Rathsman [^, on the other 
hand, do not include in their analysis the two-loop correction responsible for 
correlations within the dispersive approach, and argue that its contribution 
should be negligible. Such differences are perhaps difficult to settle when 
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working with multiparameter fits, however it seems that a more extensive 
comparison with data for different event shapes should settle this issue. 

A second important point is the fact that if the strong coupling as{Mz) is 
also included as a £t parameter, the results point to a value which is rather 
low with respect to the world average (typically as{Mz) ~ 0.110). This 
should perhaps be considered as a cause of moderate concern: it could mean 
that the current expressions for the shape function are not fully adequate to 
describe the effect of power corrections in the peak region (and thus they pull 
the £t towards lower values of Os), or it could imply that current assumptions 
about the theoretical error to be associated with as{Mz) underestimate the 
effect of higher order and power corrections. 

A hnal observation of considerable relevance for phenomenology was made 
by Salam and Wicke 1^, and it concerns the relationship between the the¬ 


oretical calculations described so far and experimental data. Briefly stated, 
the point is that while all models of power corrections are derived in mass¬ 
less QCD, event shapes are measured using the massive particles produced 
by hadronization. Salam and Wicke observed that the difference between 
massless and massive dehnitions of event shapes induces non-universal power 
corrections of the same parametric size (A/g) as conventional ones. To disen¬ 
tangle these mass effects from the universal features of soft gluon emission, for 
example when htting nonperturbative parameters, it is necessary to specify a 
scheme to connect massless QCD computations and measured event shapes. 
The procedure followed this far, which ignores mass effects in the dehnition 
of events shapes, also constitutes a possible scheme; in this scheme, however, 
different event shapes are treated differently, and non-universal, mass-related 
power corrections should be htted separately for different observables. Other 
schemes have been proposed: notably, the “T^-scheme” uses measured en¬ 
ergies Ei in the definition of the event shape, and rescales three-momenta 
by Ei/pi] in this scheme, three-momentum is not conserved, however it can 
be shown that non-universal power corrections vanish in a tube model cal¬ 
culation. Another interesting possibility is the “decay scheme”, in which all 
measured particles are forced to decay into massless particles via a Monte 
Carlo interface; not all decays are strong, nor realistic, however the method 
is interesting because it begins to address in some detail the issue of possible 
double countings, which arises when parton-level, QCD-based predictions 
for power corrections are used in conjunction with other hadronization mod¬ 
els such as those implemented in Monte Carlo evolution codes. 

Mass-related power corrections of the type described in Ref. 


23 


are 
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further logarithmically enhanced by effects related to hadron multiplicity in 
the hnal state. Model calculations using the hypothesis of local parton- 
hadron duality suggest 

A/ a\^ 4N 

6m{e) = Ce-[\og^j ; A = -^~1.6, (19) 

where 6m{e) is the non-universal power correction to the average of event 
shape e, and Cg a htted coefficient. It should be noted that conventional 
power corrections will also generically be enhanced by logarithms, however 
not much attention has so far been devoted to study their effects. 


6 Perspectives 


After more than two decades of studies, event shape distributions remain 
at the forefront of theoretical and experimental QCD analyses. In a single 
graph, a typical distribution takes us from a completely perturbative regime 
(e ~ 1), where, at least at high energy, hxed order calculations apply, to 
a region dominated by nonperturbative effects (e ~ A/g), where all power 
corrections become important. Successive improvements in our theoretical 
tools are leading to a QCD-based understanding of the entire distribution, 
including the peak region dominated by two-jet conhgurations. 

Resummed QCD amplitudes point beyond perturbation theory, and the 
effects of nonpertubative corrections can be conveniently parametrized by a 
new class of nonperturbative functions, the shape functions, which provide 
a general framework for studies of power corrections. Like parton distri¬ 
butions, they must be htted from data, however different QCD-motivated 
models suggest somewhat different functional forms |^, |^. The results 
of the dispersive approach are recovered, and the universality of the lead¬ 
ing power corrections placed in a wider perspective. Renormalon calculus 
can be merged with Sudakov resummation via dressed gluon exponentiation, 
providing a model for the shape function. Hadronization effects generate 
mass-related, log-enhanced power corrections which must be separately un¬ 
derstood: a dehning scheme for event shapes must be chosen before a detailed 
phenomenological analysis is attempted. 

Looking to the future, much work still needs to be completed to fur¬ 
ther rehne our theoretical understanding. At hxed order, studies of event 
shapes at NNLO are just beyond the horizon |^; in due course, they will be 
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supplemented by NNL resummation of Sudakov logarithms The shape 


function viewpoint could be extended in several directions: to more gen¬ 
eral event shapes in annihilation processes, which require special treatment 
within the dispersive approach, such as jet broadening 0 , energy-energy 
correlations and the D parameter |^; to other processes, such as DIS, 


where much is known already about resummed event shapes, though it was 
not covered here because of lack of space ; ultimately, to the more general 
and interesting case of QCD hard scattering. Dressed gluon exponentiation 
cannot presumably be extended beyond NLL accuracy without tackling the 
issue of double renormalon chains, in itself a question of some theoretical 
interest. In summary, many interesting questions remain open, and our un¬ 
derstanding of the interface between perturbative and nonperturbative QCD 
is likely to become deeper in the coming years. 

As for experiment, this being a brief review form a theoretical viewpoint, 
I did not discuss data (see, for example. Ref. [^). It should however be 
emphasized that one of the main reasons for theoretical interest in this area 
is the remarkable wealth and precision of the data available, over a wide 
kinematic range and for different processes. The precision of the data is 
such that they may well be able to discriminate between different current 
theoretical models of power corrections in the peak region. As the work of 
LEP collaborations appears to be winding down, it is perhaps worth stating 
once more the obvious: it is very important that data and the tools for their 
analysis should be preserved in usable form long beyond the lifetime of the 
experiments. New theoretical tools that warrant further data analysis may 
yet be developed, and good experimental work done in the past should not 
be lost. 
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